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Abstract 

This work investigates the impingement of a liquid microdroplet onto a glass substrate at 
different temperatures. A finite-element model is applied to simulate the transient fluid dynamics 
and heat transfer during the process. Results for impingement under both isothermal and non- 
isothermal conditions are presented for four liquids: isopropanol, water, dielectric fluid (FC-72) 
and eutectic tin-lead solder (63Sn-37Pb). The objective of the work is to select liquids for a 
combined numerical and experimental study involving a high resolution, laser-based interfacial 
temperature measurement to measure interfacial heat transfer during microdroplet deposition. 
Applications include spray cooling, micro-manufacturing and coating processes, and electronics 
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packaging. The initial droplet diameter and impact velocity are 80 |um and 5 m/s, respectively. 
For isothermal impact, our simulations with water and isopropanol show very good agreement 
with experiments. The magnitude and rates of spreading for all four liquids are shown and 
compared. For non-isothermal impacts, the transient drop and substrate temperatures are 
expressed in a non-dimensional way. The influence of imperfect thermal contact at the interface 
between the drop and the substrate is assessed for a realistic range of interfacial Biot numbers. 
We discuss the coupled influence of interfacial Biot numbers and hydrodynamics on the 
initiation of phase change. 

Keywords: microfluidics, drop impact, interfacial heat transfer, numerical simulation, 
thermoreflectance technique. 
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Nomenclature 

Bi Biot number [h c dokf ] 

c speed of sound [m s" 1 ] 

c p specific heat [J kg^K" 1 ] 

C dimensionless heat capacity [pc p /pic Vi i] 

d splat diameter [m] 

e distance of center of circular measurement spot from origin on r-axis [m] 

2 -1 - 1 

Fr Froude number [v do' g ] 

g gravitational acceleration [9.81 m s" ] 

2 1 

h interfacial heat transfer coefficient [Wm" K" ] 

H mean surface curvature [m" 1 ] 

H dimensionless mean surface curvature [Hdo] 

k thermal conductivity [Wnf'K" 1 ] 

K dimensionless thermal conductivity [kkf 1 ] 

M Mach number [voc -1 ] 

n number of grid points inside circular measurement spot 
p pressure [Pa] 

2 - 1 

P dimensionless pressure \pvd pf ] 

Pr Prandtl number [juc Pi ikf 1 ] 

q heat flux at the splat/substrate interface [W m" ] 

r radial coordinate [m] 

R dimensionless radial coordinate [rdo~ l ] 

Re Reynolds number [pvodojU _1 ] 
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s radius of spot [m] 

t time [s] 

T temperature [K] 

u radial velocity [m s" 1 ] 

U dimensionless radial velocity [uvd 1 ] 

v axial velocity [m s" 1 ] 

V dimensionless axial velocity [wo" 1 ] 

2 1 

We Weber number [pvo do/' ] 

z axial coordinate [m] 

Z dimensionless axial coordinate [zdo "'] 
Greek symbols 

2 1 

a thermal diffusivity [m s" ] 

P spread factor [d max do~ l ] 

8t temporal resolution available by experimental setup 

y surface energy [Jm" ] 

(|) contact angle 

|u dynamic viscosity [Pa s] 

9 dimensionless temperature [{T- min(ri j0 , 72,o)}(|^i,o-72,o|) 

p density [kg m" ] 

a stress [Pa] 

x dimensionless time [Zvo^o" 1 ] 
Subscripts 

initial 



1 drop/splat 

2 substrate 
avg average value 

c contact, interface 

i initial 

int linearly interpolated value 

/ liquid 

max maximum value 

r radial direction 

z axial direction 
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1 Introduction 

The fluid dynamics and heat transfer associated with microdroplet impingement onto a substrate 
are of considerable interest to micro-manufacturing, spray cooling, spray coating, and inkjet- 
printing [1-3]. A variety of fluids are used in such processes, including fuels in combustion, 
water and dielectric fluids for cooling, and metal droplets for rapid prototyping and electronic 
interconnects [4, 5]. 

In this work a numerical investigation of a liquid microdroplet impacting on a horizontal 
substrate at different temperatures is presented (Figure 1). The initial droplet diameter and 
impact velocity are 80 um and 5 m/s respectively and gravity is negligible. The associated 
transport phenomena are extremely complex. For instance, this problem involves fluid dynamics 
with large deformations of the droplet free surface. The simultaneous, transient heat transfer 
process involves convection in the droplet coupled with conduction in the substrate. Both the 
thermal field inside the droplet and the onset of phase change, if any, depend on the interfacial 
heat transfer coefficient between the droplet and substrate, which expresses the imperfect 
thermal contact at the interface. Our study focuses on drops of eutectic tin-lead solder (63 Sn- 
37Pb, referred as solder, hereafter), water, isopropanol and FC-72 fluorocarbon, a dielectric fluid 
used for electronics cooling. A primary objective of this work is to evaluate these liquids as 
potential candidates for a companion experimental study, currently being developed, to measure 
interfacial heat transfer coefficients and temperature history and compare the results to numerical 
simulations. 

Traditionally, numerical models targeting similar problems use strong simplifications for 
the sake of numerical tractability [3]. For instance Harlow and Shannon [6] neglected both 
viscous and surface tension effects in their modeling of a liquid droplet impacting on a flat plate. 



6 



Tsurutani et al. [7] used the simplified marker and cell method (SMAC) and employed a fixed 
grid with relatively low resolution. Increasing computing capacities have recently led to very 
convincing simulations of the impact of millimeter-size drops with the Volume-Of-Fluid method 
[8], however the ability of this technique to address micrometer-size droplet cases, where free 
surface effects are more important, is not assessed yet. Gao and Sonin [9] developed a powerful 
theoretical analysis in which order-of-magnitude approximations were made to characterize the 
associated time scales, such as the times required to remove the initial superheat, remove the 
latent heat during freezing, and subsequently cool the deposit to the ambient temperature. Two 
effects that have been shown to be significant in other studies were neglected in this formulation: 
convection effects within the droplet and thermal contact resistance at the splat-substrate 
interface [10, 11]. 

Zhao et al. [12] modeled the cooling of a liquid microdroplet, accounting for fluid 
dynamics phenomena and assuming perfect interfacial thermal contact. This group used a 
Lagrangian formulation, extending the fluid dynamics model of Fukai et al. [13] to account for 
the heat transfer process in the droplet and substrate. Wadvogel et al. [14, 15] extended this 
modeling to account for solidification and imperfect interfacial thermal contact. This modeling is 
used in this article, with the incorporation of a more stable and versatile mesh generation scheme 
Mesh2d [16], and the ability to modify the interfacial heat transfer coefficient with respect to 
time and space. 

Several studies have specifically investigated the role and importance of imperfect 
thermal contact between the substrate and the drop. This imperfect thermal contact is a critical 
parameter in the heat transfer process. Liu et al. [17] suggest that when a liquid spreads over a 
solid surface, perfect thermal contact cannot be achieved between the liquid and solid surface 
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because of the substrate surface roughness, surface tension, surface impurities, and gas 
entrapment. It is believed that heat transfer through the actual (imperfect) contact area occurs by 
conduction and, to some degree, radiation across the gas-filled gaps [3]. For molten lead 
droplets, imperfect thermal contact was experimentally observed by Bennett and Poulikakos 
[11]. Pasandideh-Fard et al. [18] and Xiong et al. [19] performed a numerical study on the 
sensitivity to contact resistance on the final diameter, overall shape and height of a solidified 
solder droplet. Their model predicted variations in solders bump height up to 20% due to 
variations of thermal contact resistance. Recently, Attinger and Poulikakos [20] compared 
experimental and numerical transient oscillations for a solidifying solder drop and were able to 
estimate the value of the interfacial heat transfer coefficient for a specific case. Although the 
investigations above have shown the importance and effects of the interfacial heat transfer 
coefficient, there is still a lack of modeling and predictive tools to determine a priori the 
interfacial heat transfer coefficient. 

This study is aimed at selecting liquids and temperature for a combined theoretical and 
experimental investigation of microdroplet fluid dynamics and heat transfer on a cooled surface. 
The laser-based technique developed by Chen et al. [21] will be modified and used to measure 
the interfacial temperature with microsecond and micrometer resolution under a spreading 
droplet. Matching the measured and calculated temperature values at the interface will allow the 
determination of the transient and local behavior of the heat transfer coefficient, which is a 
necessary step in developing predictive models for interfacial heat transfer. In this article, we 
discuss the effect of interfacial heat transfer on the heat transfer process during the impact of 
solder, water, isopropanol and FC-72 (dielectric fluid) droplets on a glass substrate. 
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2 Numerical Model 

The mathematical model is based on the Navier-Stokes and energy equations [14] applied to an 
axisymmetric geometry. All equations are expressed in a Lagrangian framework, which provides 
accurate modeling of the large deformations of the free surface and the associated Laplace 
stresses [13]. 

2.1 Fluid dynamics 

The flow inside the droplet is laminar and all thermophysical properties are assumed to be 
constant with respect to temperature. The radial and axial components of the momentum 
equation are considered along with the continuity equation. An artificial compressibility method 
is employed to transform the continuity equation into a pressure evolution equation. This method 
assumes a fluid flow that is slightly compressible, whereby the speed of sound is large, but not 
infinite. A Mach number of 0.001 is used for all simulations in this work. The derivation of the 
boundary condition at the free surface consider forces due to pressure, viscous stresses and 
surface tension [13]. The traditional no-slip boundary condition fails in the vicinity of the contact 
line because its application results in an infinite stress in the region. To circumvent this problem, 
a scheme proposed by Bach and Hassager [22] is utilized, which applies a net interfacial force 
given by the equilibrium surface tension coefficient of the joining phases. The wetting force at 
the dynamic contact line between the liquid droplet and the substrate is neglected throughout the 
analysis. The dimensional form of the fluid dynamics equations can be found in [23], with the 
expression of the stress tensor. The dimensionless equations for fluid dynamics are [14]: 
Mass conservation: 



dP 1(13 



dr M 2 yR dR y 



RU + 



dvV 
dz) y 



= 



(1) 
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where P, r, R, Z, U, V are dimensionless pressure, time, radial distance, axial distance, radial 
velocity and axial velocity, respectively. M denotes Mach number. 
Momentum conservation in radial direction: 



dU 1 d /„_ \ da R7 1 

(RCrr) — + — CT m =0 

dr RdR y RR ' dZ R 06 



(2) 



where a RR and a zz are dimensionless stress tensor terms. 
Momentum conservation in axial direction: 

(i?o- ffff ) ^ + — = (3) 

ar RdR y ' dZ Fr 

In the above equation, Fr denotes Froude number. 

2.2 Heat transfer 

The energy equation is solved in both the droplet and the substrate, according to the formulation 
in [14]. Convection and radiation heat transfer from all exposed surfaces is neglected. The 
dimensionless energy conservation equation for droplet and substrate is given by (i = 1 for 
droplet and i = 2 for substrate): 



C, 



d6 ; 



1 



dr PrRe 



R dR 



K;R 



dR 



+ - 



dZ 



K 



1 dZ 



-0 



(4) 



where C, and K { is the dimensionless heat capacity and thermal conductivity respectively. Pr and 
Re denotes Prandtl and Reynolds number respectively. 6 t is the dimensionless temperature, and 
is defined as: 



0i = 



Tj -min(r i0 ,r 20 ) 

Kl,0 _ ^2,0 



(5) 
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where T\$ and 72,0 are the initial dimensional temperature of drop and substrate respectively. 

2.3 Thermal contact resistance 

Thermal contact resistance between droplet and substrate is modeled by a thin layer of arbitrary 
thickness 8, with zero heat capacity and adjustable thermal conductivity k t [19]. The interfacial 
heat transfer coefficient can therefore be defined as h c = k/d This approach is fully compatible 
with that of Wang and Matthys [24]. The thermal contact resistance can be non-dimensionalized 
with the Biot number as [19]: 

- ~y~ ( 6 ) 

where, do is the initial diameter of the droplet. 

2.4 Initial and Boundary conditions 

The initial conditions are as follows: 

U = 0;V = -l;P = — (7) 
We y ' 

l (R,Z,O) = l;0 2 (R,Z,O) = O for solder (8) 

X (R, Z,0) = 0; 6> 2 (R, Z,0) = 1 for water, isopropanol and FC-72 (9) 

The last two initial conditions show that the solder drop is cooled upon contact with the 
substrate, while water, isopropanol and FC-72 drops are heated. 
The boundary conditions are as follows: 

dV 

^ = 0;-^- = ati? = (10) 

U = V = atZ=0 (11) 
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j j 

o RR n R + o RZ n z = -2 — n R at droplet free surface (12) 
We 





<J ZR n R + (TzzUz = -2 — n z at droplet free surface (13) 
We 

The above two boundary conditions are the balance of forces due to pressure, viscous stresses 
and surface tension at droplet free surface. 

36 30 

— l -n H -n = at droplet free surface and the substrate boundary surface (14) 

OR r dZ z 

3 Numerical Scheme 

The computational domain is discretized as a mesh of triangular elements and the numerical 
model is solved using a Galerkin finite element method. Linear shape functions are used for 
velocity and pressure. An implicit method is utilized for the integration of fluid dynamics 
equations in time, while a Crank-Nicholson scheme is used for the energy equation. Details of 
the algorithm are given in [15]. The present model uses a more robust and freely available mesh 
generator Mesh2D [16]. It has been found that Mesh2D is better than the advancing front method 
[25] in terms of the time taken to generate mesh, the allowable aspect ratio of the elements, and 
the number of elements generated. A comparison of meshes generated by two methods is shown 
in Figure 2. 

In the present work, the grid and time step independence are examined in terms of the 
position of the contact point between the z-axis and the free surface of the splat (Z c ). This study 
is carried out for a 50 |um solder droplet impacting a flat surface at 2.0 m/s under isothermal 
conditions. This corresponds to Re = 314 and We = 4 .76. The grid independence is considered 
for four increasing number of nodes in the droplet: 199, 521, 705 and 873, with a time step of 
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5x10 4 in each case. The time step independence is considered for time steps of 0.5, 1, 2, and 3 
xlO for 705 nodes in each case, with the results shown in Figure 3. As it can be seen, a time 
step of 5x1 0~ 4 and a spatial discretization of 700 nodes in the droplet are sufficient for the 
simulations. Each simulation requires approximately 6 CPU hours on a 2.4 GHz Intel-Xeon 
machine with 1 GB of RAM. 

3.1 Thermophysical properties and dimensionless numbers 

The thermophysical properties and dimensionless numbers used for simulations are given in 
Table 1 and Table 2 respectively. 

4 Results and Discussions 

Results are presented for solder, water, isopropanol and FC-72 droplets with diameter d = 80 
(j,m; velocity v = 5 m/s and values of Bi of 1, 10 or 100. This choice of Biot numbers represents 
a realistic range of values used in previous work [19]. The initial temperatures are dimensionless, 
which means that a single simulation result describes any non-isothermal impact. In case of 
solder the droplet is cooled by the substrate, so the initial dimensionless temperatures for drop 
and substrate are 1 and respectively. In the cases where water, isopropanol and FC-72 droplets 
are heated by the substrate, these corresponding values are and 1 . 

4.1 Fluid dynamics 

Figure 4 shows the spreading of a solder microdroplet with successive representations of the 
droplet shape, temperature isotherms, and streamlines. During the initial spreading stage (t < 12 
|us), the deformation of the drop is mostly influenced by inertial forces. However, in the later 
stages of spreading (t > 25 |us), inertial forces decrease and surface tension forces dominate. This 
competition between inertial and surface tension forces induces the peripheral ring visible for t = 
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20 |us, as well as a strong recoiling which results in the splashing of the solder drop. Also a 
vortex forms in the drop during recoiling (Figure 4). 

4.1.1 Comparison of spreading in all four liquids 

The temporal evolution of the spread factor /? (ratio of maximum splat diameter to initial droplet 
diameter) for all four liquids is plotted in Figure 5. The least spreading is observed with solder, 
which due to its small Weber number (Table 2) and the high viscous dissipation. The maximum 
spreading occurs with FC-72 because of its large Weber number. In general, a larger Weber 
number results in more substantial droplet spreading. 

4.1.2 Comparison with previous results 

Recently our numerical code was validated with experimental results for solder [20]. In the 
present work, numerical values of the maximum spread factor for water and isopropanol are 
compared with visualization results [26], and also with analytical expressions available in the 
literature [18]. In Table 3, we use the same parameters as in the visualization study [26]: for 
isopropanol the parameters are do = 87 um and v = 9.28 m/s (Re = 259, We = 277). For water, 
the parameters are: do = 83 |um and vo = 8.19 m/s (Re = 696, We =77). The analytical estimate 
for the maximum value of the spread factor (J3 max ) in [18] assumes that the surface energy at the 
maximum spreading equals the kinetic and surface energy before impact, less the viscous 
dissipation during impact: 



d 



max 




We + 12 



(15) 



max 



d. 







where, 0is contact angle. 
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Table 3 shows very good agreement between numerical, experimental and analytical 
results, except for the maximum spread factor of isopropanol. This may be explained by the fact 
that the lower Reynolds number related to the isopropanol impact does not fully match the 
assumption in Equation 15 that viscous dissipation is due to an established boundary layer 
between the drop and the substrate. The viscous dissipation term in Equation 15 would thus 
require a modification for Reynolds numbers lower than 500 to incorporate this effect. 

The maximum spread factor of isopropanol is greater than that for water due to its larger 
Weber number. For impact of a liquid with We > 1 , the spreading process is driven by the radial 
pressure gradient induced by the sudden velocity change at the impact location [23]. After the 
maximum spread factor is reached, the water splat recoils (Figure 5). The isopropanol case 
shows that less recoiling occurs after the maximum spread factor is attained. 

4.2 Heat transfer 

4.2.1 Effect of Biot number and drop properties 

Figure 6 shows the influence of the Biot number on the four liquids. As the Biot number 
increases, heat transfer occurs more rapidly between the substrate and the drop for all cases. This 
can be verified by the location of typical isotherms for Bi = 1 and 100. It is interesting to notice 
that the temperature gradients in the solder drop are in the radial direction, while isopropanol and 
FC-72 splats exhibit axial temperature gradients. This is due to the higher thermal diffusivity of 
the solder: during the impact, the solder drop assumes a doughnut shape, with high-temperature 
fluid continuously supplied to the center region, so that the splat periphery is rapidly cooled by 
contact with the low-temperature substrate. The association of this flow pattern and the higher 
thermal diffusivity results in radial temperature gradients. 
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In later stages of spreading for isopropanol and FC-72, the heat flux across the interface 
is mainly governed by conduction through the substrate. The low values of thermal conductivity 
(Table 1) for these two liquids result in axial thermal gradients for both. 

The occurrence of phase change (if any) can be predicted from Figure 4 and Figure 6. For 
all four liquids, when the Biot number is large (Bi = 100), the largest temperature change is seen 
at the periphery of the drop. This implies that the drop will begin to solidify or evaporate at its 
periphery first. When Bi ~ 1 the convective heat flux is reduced and the orientation of the 
isotherms in Figure 6 show that phase change will occur simultaneously over the entire contact 
surface between the droplet and substrate. 
4.2.2 Effect of droplet liquid on temperature change in splat 

The thermal diffusivities of the four liquids are listed in decreasing order in Table 1. 
Accordingly, variations of temperature inside the splat occur more rapidly for higher values of 
thermal diffusivities. This can be quantified analytically by considering the splat and substrate as 
semi-infinite bodies. The analytical solution of the transient 1-D heat conduction problem in a 
semi-infinite medium that is initially at a uniform temperature T\ s o and is put in contact at time t 
= with a semi-infinite body at another temperature T^o is [27]: 



T(z,t)-T l0 ( z \ 
— = eric — j= 



(16) 



Considering the splat as a semi-infinite medium and using the analytical approach in 
Equation 16, we can determine analytically the axial splat thickness corresponding to a 20% 
temperature change at the time corresponding to the maximum extension of spreading (Table 4). 
For Bi = 100 this thickness is 0.96, 0.66 and 0.47 um, for water, isopropanol and FC-72 
respectively. Numerically, the thickness corresponding to a 20% change in temperature at the 
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maximum extension of spreading can be determined from the simulations. These thickness 
values are 0.93, 0.64 and 0.41 |um for water, isopropanol and FC-72 respectively. In both the 
analytical and numerical approach, the time to reach the maximum spreading is obtained from 
the simulations as 24, 35 and 41 us, for water, isopropanol and FC-72 respectively. Results for 
solder are not compared because thermal gradients are in the radial direction. The comparison 
between the analytical and numerical results for the thickness corresponding to a 20% change in 
temperature gives therefore reasonably consistent results (within 15% error, Table 4), provided 
the thermal diffusivity is not too important. In the case of solder for example, thermal diffusivity 
is about 1000 times higher than isopropanol and FC-72, which induces vertical isotherms: 
therefore no comparison is possible between the analytical model and the measurement in this 
case, but only comparison between the numerical calculation and the measurement. 

5 Feasibility of experiments 

A key objective of this work is to study the feasibility of using a recent laser-based temperature 
measurement technique [21] together with the numerical simulations. This coupled study will 
provide data with unprecedented temporal and spatial resolution on the behavior of interfacial 
heat transfer during droplet impingement on a substrate. The measurement technique is a laser- 
based thermoreflectance technique that measures the temperature at the fluid-substrate interface 
[21]. This technique is being modified to probe the temperature with an improved temporal 
resolution of 1 us and a spatial resolution of 15 |um. 

The setup is shown in Figure 7. A low-power He-Ne laser and a silicon photodiode are 
used to monitor the real-time reflectivity of the interface Both the droplet and substrate have a 
temperature-dependent refractive index, with the result that temperature changes in the droplet 
and substrate induce a reflectivity change of a laser beam incident on the droplet-substrate 
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surface.. By measuring the change in intensity of light reflected from the interface, the 
temperature at the interface can be obtained. The measured temperature change AT is 
proportional to the photodiode voltage change AFand can be determined as follows [21]: 

AT = ^ *° =rAV (IV) 



V 



8R 8n l OR dn s 
dn l dT dn s dT 



where R is the reflectivity, n the refractive index, and subscripts / and s are liquid and substrate, 

3ft Oft 
respectively. Since — - is typically much less than — -, the variation of substrate reflectance is 

dT ~ ~ dT 
negligible in comparison with the variation of droplet reflectance. 

Such a non-intrusive method is an ideal candidate for local and transient interface 
temperature measurements. Matching experimental and numerical temperatures (with the Biot 
number as a parameter) will allow the determination of the interfacial heat transfer coefficient. 

5.1 Selection of liquids 

The determination of the most appropriate liquids for the measurement of heat transfer 
coefficient can be helped by simulations showing how the temperature history of the droplet- 
substrate interface evolves during impingement. Such information is shown in Figure 4 and 
Figure 6. For example, the solder splat exhibits strong variation in the interface temperature in 
the radial direction for Biot numbers in the range 100. Solder is thus a strong candidate for 
experiments studying the spatial variation of interfacial heat transfer coefficient. On the other 
hand, both isopropanol and FC-72 spread more (28% to 40%) than solder (Figure 6), which 
implies that a proportionally larger droplet area will be available for the laser measurement. It is 
worth mentioning that the spreading evolution (diameter at the interface vs. time) can be 
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measured using the same laser technique. For water, the maximum spreading of 64 \xm is 
relatively small and the radial interface temperature is moderate, which represents a combination 
of both the solder and isopropanol/FC-72 behavior. As a conclusion, isopropanol can be used for 
testing the method, while solder and FC-72 will be tested because of their practical relevance. 

5.2 Error induced by the spatial and temporal resolution of the measurement 

This numerical study also provides estimates of the needed spatial and temporal resolution of the 
laser measurement to accurately capture key features of the fluid and thermal dynamics. For 
example, the entire spreading and cooling of a solder drop takes less than 100 us with maximum 
spreading diameter of 76 |um (Figure 4). The experimental method is expected to provide an 
estimated temporal resolution of 1 us and spatial resolution of 15 um, corresponding to the 
circular laser spot at the droplet-substrate interface. It is worth estimating the error induced by 
the spatial-averaging due to the extension of the spot size, as well as the error induced by the 
time-averaging. This is shown in Figure 9, where spatially averaged temperature profiles 
(simulations of the expected outcome of the experiment) are compared with the numerically 
obtained temperatures. In the spatially-averaged profile, numerical temperatures are averaged 
within successive 15 |um spot as follows (Figure 8): 



where n is number of grid points inside the spot; 7k, int is the linearly interpolated temperature 
value at the middle of the segment joining two consecutive grid points: 7k, int = (7k + 7k+i)/2; T s \ 
and T S 2 are the temperatures at the intersection of the r-axis and the circular measurement spot 




T, 



k=\ 



(18) 




19 



(Figure 8). The area A k is determined by the intersection of the circular spot and the disk defined 
by the k t h isotherm in the r-9 plane, located at a radial distance of r k , and s is the radius of 
measurement spot. If e is distance between the spot center and the origin, A k can be expressed as 
follows [28]: 



Case I: If r k > \e — s , then 



cos 



2 2 2 \ 

i e 2 +r k 2 -s 2 ^ 



2er,, 



+ s A 



cos 



2,2 2 "\ 

A e +s -r k 



2es 



-j(-e + r k +s)(e + r k -s)(e-r k +s)(e + r k +s) 



(19) 



Case II: If r. < \e — s\ , then 



A k =xr k 



(20) 



Figure 9 compares the interface temperature obtained directly from the numerical 
simulation with the spatial-averaging procedure corresponding to a laser measurement (Equation 
18) for all four liquids on the maximum extension of their spreading. The center location of the 
spatially-averaged spots e is varied from to r c -s, with a resolution of 1 um. The actual and 
spatially-averaged temperatures are shown as solid and dashed lines, respectively. In all cases the 
spatially averaged temperatures are in very good agreement with the numerical values. The 
numerical studies provide therefore insight the uncertainty of the experiment. 

Similarly, to estimate the uncertainty induced by a measurement with a temporal 
resolution of St, corresponding to the available experimental setup, numerical temperatures are 



averaged within St as follows: 
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T avg , temporal = ^ l e T(t)dt (21) 

where, t s and 4 are the start and end time within which the temperature value is measured and t e - 

t s = St. 

Figure 10 shows the evolution of the interface temperature at R = as a function of time 
for water. The solid curve shows the numerical results while the symbols simulate a temperature 
measurement with St = 1, 5 and 10 |^s. It appears that any of these resolutions is suitable for the 
temperature measurement at times larger than 10 us after impact. However, the height of the 
very short transient temperature peak occurring in the first 3 |us after impact is underestimated by 
about 25%, even if the highest available temporal resolution available experimentally (1 |^s) is 
used. 

6 Conclusions 

A numerical investigation of the fluid mechanics and heat transfer for a liquid microdroplet 
impacting on a substrate at a different temperature has been performed. In particular the effects 
of interfacial heat transfer, droplet spreading, and temperature variation at the interface are 
assessed. The liquids investigated are eutectic lead-tin solder (63Sn-37Pb), water, isopropanol 
and FC-72. Among the liquids, the spreading of FC-72 is the largest because of its larger Weber 
number. The interfacial Biot number is shown to control the location of the onset of phase 
change: for instance phase change is shown to happen at the droplet periphery if the Biot number 
is sufficiently large (Bi > 100). The numerical results are compared with published experimental 
results as well as an elementary analytical analysis. A key objective of this work is to assess the 
feasibility of a novel laser-based measurement technique to measure interfacial temperature at 
the droplet-substrate interface, with a high temporal and spatial resolution of respectively 1 |us 
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and 15 um. To assess the feasibility of this technique, numerical results are used to predict the 
droplet spreading and temperature history. These numerical results are used to determine if the 
expected spatial and temporal limitations of the experimental technique will be sufficient to 
adequately resolve the transient temperatures at the droplet-substrate interface. The initial 
conclusions are that the experimental technique will be able to accurately capture the 
temperature history at the droplet-substrate interface, given the available temporal and spatial 
resolutions. The results also show that the eutectic solder is the best candidate to measure radial 
temperature variations, while FC-72 and isopropanol exhibit larger spreading diameters and thus 
are natural candidates for preliminary experiments. 
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9 Figure captions 

Figure 1 : Problem definition and initial mesh 

Figure 2: Comparison of meshes generated by (a) advancing front method (b) Mesh2D. 
Figure 3: (a) Grid independence study: Variation of z-axis contact point with time for different 
numbers of nodes in the splat (b) Time-step independence study: Variation of z-axis contact 
point with time for different time steps. 

Figure 4: Spreading, recoiling and splashing of solder drop. Isotherms (on left hand side) and 
streamlines (on right hand side) are shown for 0.0 to 57.6 |us (Bi = 100). Splashing occurs at this 
latter time. 

Figure 5: Evolution of spread factors with time for all four liquids. 

Figure 6: Effect of Biot number on splat shape and temperature distribution for solder, water, 
isopropanol and FC-72. Isotherms (on left hand side) and streamlines (on right hand side) are 
shown at the maximum extension of spreading of corresponding splats. 
Figure 7: A schematic diagram of proposed experimental set up. 

Figure 8: Geometry used to calculate average temperature inside a laser measurement spot in r-9 
plane. 

Figure 9: Comparison of actual and spatially-averaged temperature results simulating the 
measurement of a laser measurement spot size of 15 (J,m (Bi = 100). Solid line patterns show 
actual results while dashed line pattern denotes spatially- averaged results. 

Figure 10: Comparison of actual and time-averaged temperature results assuming a data 
acquisition time of 1, 5 and 10 |us, at the location R = on interface (for water, Bi = 100). Solid 
line pattern shows actual result while symbols denote time-averaged results. 
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10 Tables 



Table 1: Thermophysical propeties used in the simulations 



Droplet 


Density 


Thermal 


Specific 


Viscosity 


Surface 


Initial 


Thermal 




(kgm- 3 ) 


conductivity 


heat 


(Pa-s) 


energy 


dimension- 


diffusivity 






(Wm-'IC 1 ) 


(Jkg'K 1 ) 




(Jm- 2 ) 


less 
temperature 


(mV) 


Solder 


8218 


25 


238 


2.6 X 10" J 


0.507 


1.0 


i.28xio _i 


Water 


997 


0.607 


4180 


9.8 X 10" 4 


7.3X10" 2 


0.0 


1.46X10"' 


Isopropanol 


785 


0.17 


3094 


2.5 X 10" J 


2.1X10" 2 


0.0 


7.0X10" 8 


FC 72 


1680 


0.055 


1050 


6.4 X 10" 4 


1.0X10" 2 


0.0 


3.1X10" 8 


Substrate 
















F2 glass 


3618 


0.78 


557 






0.0 (solder); 
1.0 (water, 
isopropanol 
and FC-72) 
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Table 2: Dimcnsionlcss numbers for drops of different liquids 



Droplet 


Re 


We 


Pr 


Bi 


Solder 


1254.7 


32.4 


2.5 X 10" 2 


1, 10 or 100 


Water 


407.4 


27.4 


6.7 


1, 10 or 100 


Isopropanol 


128.2 


73.9 


44.6 


1, 10 or 100 


FC-72 


1050.0 


336.0 


12.2 


1, 10 or 100 



Table 3: Comparison of maximum spread factor with published results 



Droplet 


Maximum spread factor ((3 max ) 


Present work 
Numerical 


Ref. [26] 
Experimental 


Ref. [18] 
Analytical* 


Water 


2.41 


2.45 


2.44 


Isopropanol 


2.52 


2.51 


2.01 



* contact angle ()) assumed to be 90° for water and isopropanol 
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Table 4: Comparison of axial thickness for a 20% change in splat obtained by numerical 

and analytical approach 



Droplet 


Axial thickness (in urn) 
Numerical 


Axial thickness (in urn) 
Analytical 


Percentage error 


Water 


0.93 


0.96 


3.2% 


Isopropanol 


0.64 


0.66 


3.1% 


FC-72 


0.41 


0.47 


14.6% 
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11 Figures 




Figure 1 : Problem definition and initial mesh 
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Figure 2: Comparison of meshes generated by (a) advancing front method 
(b) Mesh2D. 
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Figure 3: (a) Grid independence study: Variation of z-axis contact point with time 
for different numbers of nodes in the splat (b) Time-step independence study: 
Variation of z-axis contact point with time for different time steps. 



32 




Figure 4: Spreading, recoiling and splashing of solder drop. Isotherms (on left 
hand side) and streamlines (on right hand side) are shown for 0.0 to 57.6 us (Bi = 
100). Splashing occurs at this latter time. 
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16.0 



32.0 48.0 

Time (in us) 



64.0 



80.0 



Figure 5: Evolution of spread factors with time for all four liquids. 
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19.5ns 04 

N 0.2 







Bi=100 /^~\ 


\0-4 





1 -0.8-0.6-0.4 -0.2 0.2 0.4 0.6 0.8 1 
r 



0.4 
0.3 
N 0.2 
0.1 





Water 
24.0 us 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 
r 



0.4 
0.3 
N 0.2 
0.1 





0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



Bi = 1 



Isopropanol 
35.2^s 




-1.6-1.4-1.2 -1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 O-.B-l <: -1*2 -1 -0.8-0 S-C/ -0.2 (! 2 0,1 0.6 0.8 1 1.2 1.^ 

f r 



Figure 6: Effect of Biot number on splat shape and temperature distribution for solder, water, 
isopropanol and FC-72. Isotherms (on left hand side) and streamlines (on right hand side) are shown at 
the maximum extension of spreading of corresponding splats. 
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Figure 7: A schematic diagram of proposed experimental set up. 
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n = number of grid points 
covered by spot 



• Grid points 
o Center of spot 




Figure 8: Geometry used to calculate average temperature inside a laser measurement spot 
in r-9 plane. 
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0.6 




n I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
0.2 0.4 0.6 0.8 1 1.2 1.4 

Non-dimensional radial distance 

Figure 9: Comparison of actual and spatially- averaged temperature results simulating 
the measurement of a laser measurement spot size of 15 um (Bi = 100). Solid line 
patterns show actual results while dashed line pattern denotes spatially-averaged results. 
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Figure 10: Comparison of actual and time-averaged temperature results assuming a data 
acquisition time of 1, 5 and 10 ms, at the location R = on interface (for water, Bi = 
100). Solid line pattern shows actual result while symbols denote time-averaged results. 
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